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ESTIMATES OF OPTIMAL OPERATING CONDITIONS FOR HYDRO GEN- OXYGEN 
CESIUM-SEEDED MAGNETOHYDRODYNAMIC POWER GENERATOR 
by J. Marlin Smith and Lester D. Nichols 
Lewis Research Center 

SUMMARY 

In this study the value of percent seed, oxygen to hydrogen ratio, combustion pres- 
sure, Mach number , and magnetic field strength, which maximize either the electrical 
conductivity or power density at the entrance of an MHD power generator were obtained. 
The working fluid is the combustion products of hydrogen and oxygen seeded with cesium 
hydroxide (CsOH).. Two forms of the power density are investigated: The first is that of 
the ideal theoretical segmented Faraday generator. The second is of an empirical form 
found from correlating the data of many experimenters working with generators of dif- 
ferent sizes, electrode configurations, and working fluids. 

The conductivity and power densities optimize at a seed fraction and an oxygen to fuel 
ratio in the neighborhood of 3. 5 mole percent and 7. 5, respectively. The optimum 
values of combustion pressure and Mach number depend on the operating magnetic field 
strength. 


INTRODUCTION 

In the practical design of a magnetohydrodynamic power generator for either mobile 
or base power applications, one parameter is of particular importance - size. The mag- 
net size dominates the weight for mobile applications and the cost for base power applica- 
tions. Therefore, it is desirable to operate the device with controllable conditions ad- 
justed so as to obtain the maximum power density (smallest size for a specified power 
output). 

In this study we considered a gaseous hydrogen and oxygen combustion working fluid 
seeded with cesium hydroxide (CsOH) dissolved in water. Hydrogen is desirable for 
lightweight mobile systems and was therefore chosen as a fuel because it has the highest 



energy content per unit mass. Recent discussions of a hydrogen economy (ref. 1) indi- 
cate that it may also be a candidate for large base power systems. 

The controllable parameters that are considered are the percent seed, oxygen to 
hydrogen ratio, combustion pressure, Mach number, and magnetic field strength. The 
operating values of these parameters that maximize either the electric conductivity or 
the power density at the entrance to the generator are presented. 

Two forms of the power density are considered. The first is that of the ideal theo- 
retical segmented Faraday generator. However, inasmuch as the presently measured 
power density of MHD generators is considerably below this ideal value, particularly in 
the region of large Hall parameters, an empirically derived formula was also considered. 
This empirical formula was found by correlating the data of many experimenters working 
with different generator sizes, electrode configurations, and working fluids. This cor- 
relation indicated that the broad base of data could be fitted by the ideal power density 
times a correction factor which depends solely on the Hall parameter. 

A previous optimization study was reported in reference 2, where power density at 
the entrance was optimized as the first step in optimizing the entire MHD channel. In 
that case, however, an idealized hydrocarbon fuel burned in air was considered, the 
combustion gases were seeded with potassium sulphate, and only stoichiometric com- 
bustion was considered. 


SYMBOLS 


A conducting cross sectional area 
B magnetic field strength 
E electric field 
e unit electric charge 
h distance between Faraday electrodes 

I electric current 

j electric current density 

K load parameter 
k Boltzman constant 
l distance between Hall electrodes 
M Mach number 
m mass 

n number density 
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O/F oxygen to hydrogen weight ratio, eq. (16) 

P power density 

p gas pressure 

(Q e i) momentum transfer cross section between electron and 1 species 
R electrical resistance 

T gas temperature 

Tj temperature corresponding to ionization potential 
u gas flow velocity 

V voltage 

fi Hall parameter 

v total electron collision frequency for momentum transfer 

ii. 

*ei momentum transfer collision frequency between electron and i cn species 

a electrical conductivity 

cp tan a, where a is defined on fig. 2 

Subscripts: 

crit critical 

D digital conducting wall 

e electron 

eff effective 

F Faraday (used in appendix) 

H Hall (used in appendix) 

IV 

i i in species 

id ideal 

m maximum 

n neutral particle 

S shorting 

s seed particle 

x Hall 

y Faraday 


ANALYSIS 


Performance Parameters to be Optimized 

The purpose of this study is to establish the operating conditions for a hydrogen- 
oxygen, cesium seeded MHD generator based on the optimization of the generator inlet 
conditions. The optimization of three quantities is considered. The first quantity is the 
electrical conductivity since its optimization will give the maximum current density for a 
given electric field. 

The second quantity is the ideal power density. For the Faraday type of generator 
(ref. 3) that will be considered here, the ideal power density (obtained when the load 
resistance matches the generator internal resistance) is 

P id = J ctu2b2 (1) 

where a, the electrical conductivity, is 


m e v 

and where n is the electron number density, e is the unit electric charge, m is the 
electron mass, v is the electron collision frequency for momentum transfer, u is the 
gas flow velocity, and B is the magnetic field strength. 

Because a and u are basically independent of B, equation (1) indicates that the op- 
timum power density is obtained at maximum B. However, one finds that in practice 
the generator performance rapidly deteriorates as the Hall parameter /3 exceeds ap- 
proximately unity, where 


0 = 


eB 

m e v 


( 3 ) 


It is shown in the next section that a reasonable measure of the performance dete- 
rioration with /3 is given by 


P H = I cru 2 B 2 -- + L 
eff 4 P 

4 i + /r 


0 > l 


( 4 ) 


P eff ~ P id ** 1 
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where P gff is the effective power density. 

This is the third quantity which will be optimized. 


Calculation of Parameters 

Effective power density . - Ohm’s law, appropriate for an MHD working fluid in the 
presence of a magnetic field, is, in vector form (ref. 3, eq. 10. 56) and for the orien- 
tation shown in figure 1, 

j = cr(E + uxB)-/3j xb (5) 

Equation (5) can be written in component form: 

J y + fl x -a(iiB + Ey) (6a) 

j x - /3j y = ctE x (6b) 

where j is the Hall current density, j is the Faraday current density, E is the Hall 
x y x 

electric field, and Ey is the Faraday electric field. The power density generated by 
such a working fluid is 


P = j x (-E x ) + jy(-Ey) (7) 

MHD generators are built with electrodes segmented in the x direction, and three types 

of connections are commonly considered (see fig. 2): (1) Faraday, where j = 0, and 

for which E__ = -/3(uB + E„); (2) diagonal, where E„ = <pE„ with <p = tan a where a 
a y y x 

is the angle the wall segment makes with u x B; and (3) Hall, where Ey = 0. 

For each mode equation (7) shows that the power density depends on the value of E x 

that it is possible to develop. Equation (6b) shows that the value of E x for any j is 

proportional to jy. Since electrode voltage drops limit the current j y , they will also 

limit E . For a given j any current j_ will short out E and limit its value, 
x y x x 

Hence, E x is a measure of the generator’s ability on the one hand to provide j and on 
the other hand to prevent j x . The quantity 


k h M = 


- E xl^ 

uB <p + i 3 


( 8 ) 


where E x is evaluated at the condition E y = <pE x can be used as a measure of gener- 
ator performance because it is the fraction of the ideal axial electric field which a 
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generator attains. As shown in the appendix the power densities of the three types of 
generators (when loaded so as to give maximum power) are 


1 + K H (O)0 2 

2 — Faraday (9) 

i^ 2 

2 

— = K w (<p) — Diagonal conducting wall (10) 

i«v (X + A i + p 2 ) 

4 


k h (o) 


1+0 


Hall 


( 11 ) 


these results are a generalization of those given in reference 4. These forms allow for 
internal current leakage and electrode voltage drops by evaluating Kg from experimen- 
tal results. 

The results for Kg from many different experiments are shown in table I. The 
data are taken from generators with different working fluids as well as different elec- 
trode connections. The references from which the data were obtained and, in some 
cases, the data used in determining Kg are also shown. 

The data corresponding to the highest values of Kg are shown in figure 3. Also 
shown in the empirical relationship 


K „(<p) = - 0 > 1 

H 0 

= 1 0 < 1 


(12) 


which correlates the data with reasonable accuracy. 

This correlation agrees with a result obtained in reference 5 for inert gas gener- 
ators when fluctuations are taken into account. Thus substituting equation (12) for |3 > 1 
into equation (9) gives 



1+0 


1+0 


u 2 B 2 
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Then defining 


results in 


eff 


= a 


1 +/3 

1 + J8 2 


a eff 

lim/3—°° 


a 


1 

P 


which is the expression derived in reference 5 (see eq. (19) and fig. 6) for = 1. 

Although fluctuations of the type discussed in reference 5 are only predominate in inert- 
gas, nonequilibrium, MHD generators, it appears that a similar degradation of perfor- 
mance can occur in combustion MHD generators. 

Collision frequen cy for momentu m transfer . - To evaluate the conductivity and Hall 
parameter (eqs. (2) and (3)), we need the electron collision frequency for momentum 
transfer. This quantity is defined in reference 3 as 




(13) 


where the sum extends over all species except electrons and where n. is the number 
density of the i** 1 species, (Q o{ > is the momentum transfer cross section between an 
electron and the r species, k is the Boltzmann constant, and T is the gas temper - 
ature. The quantity (Q e j> is calculated from the monoenergetic beam cross section, 
Q ei (by ref. 6) as 



where v is the relative velocity between the collision partners. One of the greatest 
sources of uncertainty in combustion MHD generator analysis is in the lack of good cross 
section data in the electron energy range of interest (2000 to 4000 K). 

Before proceeding to a discussion of the cross section data used in this report, it is 
interesting to note the mole fractions of various species in the hydrogen -oxygen -cesium 
chemistry. These are listed in table H for conditions typical of those analyzed in this 
report. A discussion of the chemical analysis is presented in the next section. The im- 
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portant thing to note is that the mole fraction of water (HgO) is nearly an order of mag- 
nitude greater than any of the other species. Furthermore, its cross section is larger 
than for the other species (as shown in fig. 4), except for cesium which only appears in 
trace amounts. The result is that the collision frequency is grossly dominated by H^O. 
Thus, the only cross section that must be established with some care is that of HgO- 

The value for the monoenergetic cross section for HgO used in this report was taken 
from page 9 of reference 7. The momentum transfer cross section obtained by numer- 
ical integration of equation (14) is plotted in figure 5 along with the data used by five 
other researchers in the MHD field. It is seen that there is a wide variation in the 
values used. The value used in this report is the largest; hence, the results obtained 
from the optimization of electrical conductivity and power densities are on the conserv- 
ative side. 

Sources of other cross sections used in this report are listed here: 

(1) Cesium - Cesium momentum transfer cross section was previously evaluated in 
reference 6. 

(2) Atomic hydrogen (ref. 8, p. 132) - The curve used was an extrapolation of the 
results of Temkin and Lamkin. 

(3) Hydroxyl ion (OH) - A readily available monoenergetic cross section was not 
found. Therefore, we arbitrarily took a value between the two points given in refer- 
ence 9 and made the cross-section constant. 

(4) Molecular hydrogen - The momentum cross section is given directly in refer- 
ence 8 (p. 127). 

(5) Molecular oxygen - The experimental results given in reference 8 (p. 144) were 
used. 

(6) Atomic oxygen - The curve labeled beam experiment, Bates and Massey, and 
Lin and Kivel, was used (ref. 8, p. 133). 

(7) Ions - The collision frequency given by equation (29) of reference 6 was used. 


Method of Solution 

The equilibrium chemistry of hydrogen-oxygen-cesium combustion was calculated 
using the rocket option of reference 10. This program calculates the mole fractions of 
the 31 species listed in table H as functions of combustion pressure and isentropic ex- 
pansion ratio. Because collision cross sections were not known for all of the species, 
the unknown were assigned the cross sections given here: 

(1) Cesium cross section - Cs, CsO, CsOH, CSg, CSgO, CSgOgHg 

(2) HgO cross section - HgO, HOg, HgOg 
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<g 

(3) Hg cross section - Hg 

(4) Og cross section - Og, Og 

(5) H cross section - H 

(6) O cross section - O 

(7) Ion cross section -Cs + , CsOH + , H + , H", 0 + , O", OH + , OH“, 0 2 " 

(8) OH cross section - OH 

On this basis the conductivity and power densities given by equations (1), (2), and 
(4) are readily calculable as functions of the following variable parameters: 

(1) Seed fraction - 1 part CsOH dissolved in two parts water. The CsOH is intro- 
duced into the combustion chamber as a solid dissolved in water (as a liquid). Both were 
introduced at room temperature. The percent seed is defined as 

Moles (CsOH + 2H«,0) 

% seed = - — 100 (15) 

total moles 

(2) Oxygen to hydrogen weight ratio - 

weight of 0 0 

O/F = i (16) 

weight of Hg 

The oxygen and hydrogen were introduced into the combustion chamber as room temper- 
ature gases. 

(3) Combustion chamber pressure - The assigned value of combustion chamber 
pressure. 

(4) Area ratio - The ratio of inlet area of the MHD duct to the throat area of 
combustor. 

(5) Magnetic field - The assigned value of magnetic field strength. 

The question investigated in this study is what values of these five parameters should 
be chosen to optimize the conductivity or the power densities. The range of values 
studied are seed fraction all values; O/F, 4 to 12; combustion pressure, 0. 5 to 4 mega- 
pascals; area ratio, 00 subsonic to 10 supersonic; and magnetic field, 0 to 10 teslas. 

The procedure adopted was to investigate the effect of varying one parameter at a 
time to find an operating value that optimizes the conductivity or power densities over 
the range of values of the other parameters. This procedure is successful primarily be- 
cause the optimum values vary slowly as a function of seed fraction, O/F, and Mach 
number. 
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RESULTS AND DISCUSSION 


Optimum Seed Fraction 

The influence of seed fraction on the optimum values of conductivity and power den- 
sities is evaluated by plotting conductivity and power densities as a function of seed frac- 
tion for various values of one of the other four quantities (O/F, combustion pressure, 
area ratio, magnetic field) while holding the other three values at an arbitrarily chosen 
value which is near the midpoint of their range. These fixed values were O/F = 8 
(stoichiometric); combustion pressure, 1. 03 megapascals, and area ratio, 1.30 (at 

M = ~1.6). Because conductivity is independent of magnetic field and the ideal power 
o 

density varies as B (eq. (1)), there is no optimum value as a function of magnetic field 
for these two parameters. Hence, the influence of seed fraction on magnetic field varia- 
tion involves only the optimization of the effective power density. This effect differs for 
low and high values of magnetic field. For low magnetic field (i. e. , /3 < 1) the effective 
power density (eq. (4)) equals the ideal power density (eq. (1)). Hence, the influence of 
seed fraction on effective power density at low magnetic field is the same as it is on ideal 
power density. Therefore, the fixed value of B for all computations of seed effect was 
set at the upper range of magnetic field, namely, 10 teslas. 

Figure 6 shows the results. The only important thing to note is that, except for the 
curves at the very low value of O/F = 4 shown in figure 6(a) (which these curves show 
is far below the optimum O/F ratio), there is a slow variation of all three parameters 
as a function of percent seed. As a result, at approximately 3. 5 mole percent seed, all 
the curves are within a few percentage points of their optimum value. Therefore, a seed 
fraction of 3. 5 percent is chosen as the optimal seed fraction and is used to determine 
the optimal values of the other parameters. 

The power densities optimize in these figures primarily because the conductivity op- 
timizes with increasing seed fraction. This can be readily illustrated for an inert gas 
containing an ionizable seed. The electron number density n g is related to the seed 

number density n by Saha's law (ref. 3). 
s 


n e “ exp 



(17) 


where T- is the temperature corresponding to the ionization potential. Ignoring for 
simplicity electron- ion collisions (weakly ionized gas), we obtain a collision frequency 
(given by eq. (13)) that is proportional to 


v 01 n n ( Q en ) +n s (Q es ) 


(18) 
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Substituting equations (17) and (18) into equation (2) shows that the conductivity is pro- 
portional to 


a cc 



n n<^n> + n s< C W 


(19) 


At low seed fractions (n n ( Q^ n ) » n g ( Q gg ) ) the conductivity increases as the square 
root of the seed fraction. The equation reaches a maximum when n fl ( Q gn ) = n s (Q es > . 

A more detailed study of this effect in inert gas systems is contained in reference 11. 

In the combustion MHD system considered here , the location of the maximum is not 
as straightforward since it is masked by the chemistry. In addition the alkali metal seed 
is introduced in the form of CsOH dissolved in water, and the water acts as a diluent 
which decreases the gas temperature. From equation (19) this obviously results in a de- 
crease in conductivity and hence also in power density. 


Optimum Oxygen -Hydrogen Ratio 

Having set the optimum seed fraction at 3. 5 percent we now proceed to optimize the 
O/F. The magnetic field remains fixed at 10 tesla as discussed in the previous section. 
The two remaining variables are the combustion pressure and area ratio. In figure 7 
conductivity and the power densities are plotted as functions of O/F with area ratio as 
a parameter for three combustion pressures. In these nine figures all the curves are 
within 10 percent (and most are within a few percent) of their maximum at an O/F of 
approximately 7. 5. Therefore, an O/F of 7. 5 is chosen as the optimal O/F and is 
used to determine the optimal values of the other parameters. 

As before, the power densities show an optimum in these figures primarily due to 
the influence of O/F on the conductivity. On either the fuel rich or oxygen rich side of 
stoichiometric (O/F = 8) there is a surplus of particles that have nothing with which to 
react and hence act as a diluent and decrease the temperature and hence the conductivity 
(eq. (19)). However, the conductivity plots of figure 7 indicate that the peak in the con- 
ductivity occurs on the fuel rich side of stoichiometric. This is due to the chemistry of 
the system and occurs because cesium is a strong reducing agent. Therefore, as O/F 
increases an increasingly greater amount of CsOH and CsO is formed. This results in 
less atomic cesium, which produces the ionization, and thus results in a lower free elec- 
tron density and hence in a lower conductivity. 
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Optimum Mach Number 


With the optimum seed fraction set at 3. 5 percent and O/F at 7. 5, we now investi- 
gate the influence of area ratio (or equivalently the influence of Mach number) on the op- 
timum values of conductivity and power densities. From the conductivity plot of fig- 
ures 6(c) and 7, it is obvious that the conductivity decreases with increasing Mach num- 
ber and is therefore a maximum in the combustion chamber (where M = 0). This is due 
to the strong temperature dependence of the conductivity. 

In figure 8 au (which is proportional to ideal power density) is plotted versus Mach 
number for various combustion pressures. The maximum occurs at a Mach number of 
approximately 2. 5. The Mach number for maximum power density shows only a slight 
decrease with increasing combustion pressure. The maximum occurs as a result of the 
interplay between the rapidly decreasing value of conductivity and the increasing value of 
velocity with Mach number. 

Also in figure 8 the effective power density P^ is plotted versus Mach number for 
various combustion pressures. The maximum values of power density for the various 
pressures form a locus such that a maximum value of effective power density occurs for 
a chamber pressure between 1 and 20 megapascals. The effect of pressure is discussed 
in greater detail in the next two sections. 


Optimum Combustion Pressure 

The conductivity plots of figures 6(b) and 7 show that the maximum conductivity oc- 
curs at the lowest combustion pressure. This results from the fact that the temperature 
is only slightly affected by combustion pressure for the pressure range considered here- 
in so that for the same seed fraction equation (19) reduces to 


cr <x 



( 20 ) 


that is, the conductivity decreases approximately as the inverse square root of the pres- 

p 

sure. The oru plots of figures 6(b) and 7 indicate that the ideal power density follows a 
similar trend, due to the fact that the velocity is only slightly affected by combustion 
pressure. Hence the ideal power density depends on combustion pressure only through 
the conductivity. 

For high combustion pressure (/3 < 1) the effective power density reduces to the ideal 
power density. The power density decreases with increasing combustion pressure as 
does the conductivity. For the constant area ratio case of figure 6(a) (top plot) the 
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temperature is only slightly affected by combustion pressure. Therefore, from equa- 
tions (3), (13), and (18) the Hall parameter is approximately proportional to 

^ccB (2i) 

P 

Hence as /3-*°o from equations (4), (20), and (21) yield 

P eff 03 (22) 

and for /3 » 1 , the effective power density increases with increasing combustion pres- 
sure. Therefore, the effective power density goes through a maximum with increasing 
combustion pressure for a given magnetic field strength (shown in figs. 6(b) (bottom plot) 
and 8). Figure 8 also shows that the Mach number for optimum effective power density 
decreases with decreasing combustion pressure. Since for a given combustion pressure 
the static pressure in the MHD channel is determined by the Mach number, this result 
indicates that there probably is an optimum static pressure for a given magnetic field 
strength. Or conversely, for a given magnetic field strength there is an entrance static 
pressure and hence combustion pressure and Mach number which optimizes the effective 
power density. 


Maximum Effective Power Density and Its Consequences 

As noted from figure 8 the maximum values of P for various combustion pres- 
sures form a locus of points such that a maximum value occurs as a function of combus- 
tion pressure. This is a consequence of the empirical correlation of figure 3. For the 
magnetic field strength of 10 teslas used, this maximum occurs between combustion 
pressures of 1 and 2 megapascals. It also may be noted that this maximum value then 
sets the Mach number at which one should operate. It is therefore concluded that mag- 
nitude of the magnetic field strength sets the combustion pressure and Mach number at 
the entrance of the MHD generator for maximum effective power density. 

The optimum combustion pressure and Mach number for maximum effective power 
density are plotted versus magnetic field strength in figures 9 and 10. The optimum 
combustion pressure increases nearly in direct proportion to the magnetic field strength, 
and the Mach number is nearly constant at a value of 2. The static pressure correspond- 
ing to these operating combustion pressures and Mach numbers is plotted in figure 11 as 
a function of magnetic field strength. This static pressure is nearly proportional to the 
magnetic field strength which indicates that the maximum effective power density should 
occur at nearly constant /3 (see eq. (21)). 
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This observation is confirmed by the plot shown in figure 12, wherein /3 at the op- 
erating conditions of figures 10 and 11 is plotted versus magnetic field. That /3 for 
maximum power density is constant at a value of 1 should not be totally unexpected. The 
gas temperature and velocity for a given expansion Mach number are relatively insen- 
sitive to combustion pressure. Since the optimum Mach number is relatively constant 
with operating conditions, that is, magnetic field, it follows that the gas temperature and 
velocity are relatively insensitive to operating conditions. Therefore, for a given value 
of magnetic field the conductivity from equation (20) is proportional to 

a cc _L 

v? 

and the Hall parameter from equation (2 1) is proportional to 

0 cc I (23) 

P 

Therefore, 

a o c V/3 (24) 

Substituting equations (23) and (24) into equation (4), we find that the effective power den- 
sity is proportional to 


P eff cc V/3 (25) 

1 + /3 2 

Maximizing equation (25) with respect to /3, we find that /3 must satisfy 

1+1/3 -1/3 2 -I/3 3 = -(/3 - 1)(— Z3 2 + 2/3 + 1^ = 0 (26) 

2 2 2 2 \2 2 / 

Equation (21) is satisfied when /3 = 1. (The other two roots are negative. ) The value 
/3 = 1 is a function of the empirical model chosen (see eq. (12)). If a different expres- 
sion had been chosen, /3 would take on a different optimum value. 

Additional calculations were made to evaluate the variation of the maximum effective 
power density with magnetic field strengths. The results are shown in figure 13. The 
effective power density increases monotonically with magnetic field, as does the ideal 
power density. The corresponding open circuit voltage uB is plotted versus magnetic 
field in figure 14. 
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CONCLUDING REMARKS 


In this study we have determined the values of percent seed, oxygen to fuel ratio, 
combustion pressure, Mach number, and magnetic field which optimize the values of 
either conductivity or power densities at the entrance of an MHD generator using a 
cesium seed hydrogen- oxygen working fluid. Although the utilization of these values as 
operating parameters for an MHD generator should result in the shortest length device, 
it is not clear how these results would be affected by percent power extraction and power 
level of the device. These questions are currently being investigated. 


SUMMARY OF RESULTS 

The determinations of near optimum operating conditions for a hydrogen -oxygen, 
cesium-seeded MHD power generator have been shown based on the maximization of 
either electrical conductivity or power density at the entrance of the generator. Two 
forms of the power density were investigated. The first was the ideal value derived 
from theory and the second was an effective value based on an empirical formula obtained 
from the analysis of a broad range of experimental data. 

The conductivity and power densities maximize within a narrow range of molar seed 
fractions and oxygen to fuel weight ratios. However , within the range of conditions in- 
vestigated (percent seed 0 to 10; O/F = 4 to 12; combustion pressure, 0. 517 to 4. 14 
mPa; Mach 0. 5 to 3) the conductivity and power densities are slowly varying with per- 
cent seed and oxygen to fuel ratio so that the near maximum values of conductivity and 
power densities occurred at 3. 5 percent seed and an oxygen to fuel ratio of 7. 5. There- 
fore, these values appear to be the near optimum operating conditions for percent seed 
and O/F ratio. 

Using these values of percent seed and oxygen to fuel ratio the effect of Mach num- 
ber, combustion pressure, and magnetic field strength upon the optimum values of con- 
ductivity and power densities was investigated. The maximum electrical conductivity 
was found when the gas is at the combustion temperature or at M = 0. This reflects the 
strong temperature, that is, free electron density, dependence of the conductivity. Since 
the temperature is highest at M = 0 so is the conductivity. The ideal power density 
maximum occurred at M ^ 2. 5 nearly independent of combustion pressure. Since the 
ideal power density P- d is proportional to cru^, the maximum occurs as a trade-off 
between the rapid increase in velocity with Mach number and the decrease of conductivity 
due to decreasing temperature with Mach number. Both the conductivity and ideal power 
density were monotonic decreasing functions of increasing combustion pressure. In both 
cases this is due to the nearly inverse, square root dependence of conductivity on pres- 
sure. For the transport theory formulation of the conductivity used in this report, the 
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conductivity is independent of magnetic field while the ideal power density increases with 
the square of the magnetic field. 

The optimum values of the effective power occurred at strongly coupled values of 
Mach number, combustion pressure, and magnetic field strength. These quantities op- 
timized the effective power density such that the Hall parameter was unity. This result 
is sensitive to the empirical correlation chosen. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 18, 1976, 

506-25. 



APPENDIX - EFFECTIVE POWER DENSITY FOR 
VARIOUS ELECTRODE CONFIGURATIONS 

MHD generators can be viewed as two coupled generators with equivalent circuits 
for Faraday and Hall connections: 



where h is the distance between Faraday electrodes and l is the distance between Hall 
electrodes. The conducting cross section area is A, resistances R, the voltage V, 
and the current I. The subscripts correspond to Hall H and Faraday F generators. 
These equivalent circuits can, in turn, be replaced by a second pair of equivalent 
circuits : 
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These coupled voltage equations for each generator are 


K mF h (>* + V y) ' ^ + iVp) = V F 

K mH^ (“B - y) - + A V = V H 


(Al) 


(A2) 


The load voltages for the generators can be determined by solving these equations: 



It is now possible to derive expressions for the power density for two types of 
generators. 

First, consider a generator with the condition — °o. This is called a Faraday 
generator. In this case there is no current flowing in the Hall load resistance, and the 
Faraday voltage is 



huB 



The maximum power occurs when the load resistance R p is 


Rjr - 


K mF R iF 

1 + 


(A 5) 


(A 6) 
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For this loading of the generator, the generator power density is 


P = 


louV 


av f +k (k AV h 

K m F ( 1+K mF K mH^ 1 + ^ 



2 


(A7) 


Second, consider a generator with the condition Rp = 0. This is called a Hall gen- 
erator. Since there is no voltage across the Faraday load resistance, the power is gen- 
erated entirely in the Hall circuit, and the Hall voltage is 


K 


AV 


H 


'H 


mH 


Zuy3B 


ZuB/3 


1 + K. 


mH 


r 


iH' 


R 


H , 


(A8) 


The maximum power occurs when the load resistance Rjj is 

Rr = K mH R iH ( A9) 

The power density for this generator loading is 


P„ = I oru 2 B 2 
H 4 



(A10) 


The diagonal conducting wall generator can be analyzed by using one equivalent 
circuit: 
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which, in turn, can be written as 



where 


K 


mD 


R SP 

R SD + R iD 


R _ (i + ni 

K iD 


CT (1 + (p )A d 


The load voltage is 


"mD^ 


V D = 


a + 


| ZuB - AV d 


1 + K, 


mD 


/R ip' 
R t 


The maximum power occurs when the load resistance is = K ^^ R.^, 
density becomes 


P = I ctu 2 B 2 
D 4 


AV 


K 


D 


mD 


ZuB ( 

k ±A 


Up + &Y 


K. 


mD 


(1 + /)(1 + iS 2 ) 
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(All) 


and the power 


(A 12) 



If electrode drops are negligible, then the power density divided by the ideal power 
density is 


K. 


mF 


1 + W 

I** 2 


Faraday 


I au 2 /3 2 
4 


= < 


V 

1 + Z3 2 


K. 


mD 


Up ± ff 2 ) 




(l + fyl + f) 


Hall 


Diagonal conducting wall 


Even though K m jj and are defined in terms of resistances, they can be eval- 

uated from the voltage measurements (eq. (A3)). 


'H 


= K 


Ivl/3B 


mH 


when Rp = 0 and for the diagonal conducting wall (eq. (All)) 

V D (1 + cp 2 ) _ ^ 

ZuB( (p + fi) mD 


Thus, if 


!%(<?) =ilEZl(l± 


uB 


<p+ P 


where it is understood that E y = <pE x and K mF is set equal to 1, although all of the 
power densities can be expressed in terms of K H (<p) as 
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Faraday 


r \ + K h (0 )(3 2 


1 + 0 * 


< 


1 + $ 


Hall 


k h (<p)(<p+P ) 2 

— Diagonal conducting wall 

J1 +F)(1 + cp 2 ) 



REFERENCES 


1. Gregory, Derek P. : The Hydrogen Economy. Scien, Amer. , vol. 228, no. 1, 

Jan. 1973, pp. 13-21. 

2. Carter, C. ; and Heywood, J. B. : Optimization Studies on Open -Cycle MHD 

Generators. Proc. IEEE, vol. 56, no. 9, Sept. 1968, pp. 1409-1419. 

3. Sutton, George W. ; and Sherman, Arthur: Engineering Magnetohydrodynamics. 

McGraw-Hill Book Co. , Inc. , 1965. 

4. Celinski, Z. N. : Electrical Equivalent Circuits of DC MHD Generators. Proc. of 

the Symp. on Magnetohydrodynamic Electrical Power Generation, vol. 1, Int. 
Atomic Energy Agency 1966, pp. 323-331. 

5. Solbes, A. : Instabilities in Non-Equilibrium MHD Plasmas, A. Review. AIAA 

Paper 70-40, Jan. 1970. 

6. Heighway, John E. ; and Nichols, Lester D. : Brayton Cycle Magnetohydrodynamic 

Power Generation with Nonequilibrium Conductivity. NASA TN D-2651, 1965. 

7. Brown, Sanborn C. : Basic Data of Plasma Physics. John Wiley & Sons, Inc. , 1959. 

8. McDaniel, Earl W. : CoHision Phenomena in Ionized Gases. John Wiley & Sons. , 

Inc. , 1964. 

9. Mori, Y. , et al. : Experimental Study on Electrical Conductivity of Combustion Gas 

Plasma and Phenomena Near Electrode. Fifth Int. Conf. on Magnetohydrodynamic 
Power Generation, vol. I, Int. Atomic Energy Agency, 1971, pp. 11-26. 

10. Gordon, Sanford; and McBride, Bonnie J. : Computer Program for Calculation of 

Complex Chemical Equilibrium Compositions, Rocket Performance, Incident and 
Reflected Shocks, and Chapman -Jonguet Detonations. NASA SP -273, 1971. 

11. Smith, J. Marlin; and Shair, F. H. : Theoretical Optimum Seed Concentrations in 

Slightly Ionized Nonequilibrium Plasmas. Phys. Fluids, vol. 9, no. 8, Aug. 1956, 
pp. 1615-1617. 

12. Wu, Y. C. L. ; et al. : Two -Terminal Connected Open Cycle MHD Generators. 

Ninth Symposium on Engineering Aspects of Magnetohydrodynamics. Univ. of 
Tennessee Space Inst. , 1968, pp. 85-98. 

13. DeMontardy, A. ; and Pericart, J. : Electrical Performance of a Duct with Seg- 

mented Electrodes Under Various Conditions. Proc. IEEE, vol. 56, no. 9, Sept. 
1968, pp. 1547-1555. 

14. Teno, J. ; et al. : Hall Configuration MHD Generator Studies. Proc. of the Symp. 

on Electricity from MHD. vol. HI, Int. Atomic Energy Agency, 1966, pp. 603- 
614. 


23 



15. Sonju, O. K. ; Teno, J. ; and Brogan, T. R. : Comparison of Experimental and 

Analytical Results for a 20-MW, Combustion -Driven, Hall Configuration MHD 
Generator. Eleventh Symp. on Engineering Aspects of Magnetohydrodynamics. 

D. G. Elliott, ed. , California Institute of Technology, 1970, pp. 5-10. 

16. Dicks, J. B. ; et al. : Continuation of Diagonal Conducting Wall Generator 

Research. Tenth Symp. on Engineering Aspects of Magnetohydrodynamics. 

M.I.T, 1969, pp. 183-193. 

17. Wu, Y. C. L. ; et al. : Theoretical and Experimental Studies of MHD Power 

Generation with Char. Twelfth Symp. on Engineering Aspects of Magnetohydro- 
dynamics. Argonne National Laboratory, 1972, pp. H. 1. 1 - H. 1. 14. 

18. Holt, J. F. ; et al. : Experiments with the KTVA-1 Open Cycle MHD Generator 

System. Thirteenth Symp. on the Engineering Aspects of Magnetohydrodynamics. 
Stanford Univ. , 1973, pp. n. 2. l-II. 2. 7. 

19. Mori, Y. ; et al. : Experiments on MHD Generation with ETL Mark II. Proc. of 

the Symp. on Electricity from MHD. vol. IV, Int. Atomic Energy Agency 1968, 
pp. 2761-2768. 

20. Mattsson, A. C. J. ; et al. : Energetics 6: Magnetohydrodynamics Power. Mech. 

Eng. , vol. 88, no. 11, Nov. 1966, pp. 38-41. 

21. Louis, J. F. ; Lothrop, J. ; and Brogan, T. R . : Fluid Dynamic Studies with a Mag- 

netohydrodynamic Generator. Phys. Fluids, vol. 7, no. 3, Mar. 1964, pp. 362- 
374. 

22. Bertolini, E. ; et al. : Closed Cycle MHD Subsonic Power Generation Experiments at 

Frascati. Thirteenth Symposium on the Engineering Aspects of Magnetohydro- 
dynamics. Stanford Univ. , 1973, pp. I. 1. 1-1. 1. 11. 

23. Miyata, M. : Effects of JXB Force on the Shock Tube Driven Nonequilibrium Faraday 

Generator. Thirteenth Symp. on the Engineering Aspects of Magnetohydro- 
dynamics. Stanford Univ. , 1973, pp. I. 5. 1-1. 5. 6. 

24. Shioda, S. ; et al. : Experiments in a Faraday MHD Generator with a Supersonic 

Seeded Helium Flow within a Relaxation Region. Twelfth Symp. on the Engineer- 
ing Aspects of Magnetohydrodynamics. Argonne National Laboratory, 1972, 
pp. 1.7. 1-1.7. 10. 

25. Bertoline, E. ; et al. : Closed Cycle MHD Experiments with a Large Blow Down 

Facility. Ninth Symp. on the Engineering Aspects of Magnetohydrodynamics. 

Univ. of Term. Space Institute 1968, pp. 130-141. 


24 



26. Decher, Reiner; and Kerrebrock, Jack L. : Electrode Wall - End Loop Shorting in a 

Non- equilibrium MHD Generator. 9th Symp. on the Engineering Aspects of Magne- 
tohydrodynamics. Univ. of Tennessee Space Institute, 1968, pp. 142-155. 

27. Zauderer, Bert. ; and Tate, Eric. : Electrical Characteristics of a Linear, Non- 

Equilibrium, MHD Generator. Proc. IEEE, vol. 56, no. 9, Sept. 1968, 
pp. 1535-1547. 

28. Zauderer, B. ; and Tate, E. : Experiments in a Large Non -Equilibrium MHD 

Generator with Cesium Seeded Noble Gases and Heated Electrodes. Fifth Int. 

Conf. on Magnetohydrodynamic Power Generation, vol. II, Int. Atomic Energy 
Agency, 1971, pp. 431-446. 

29. Buznikov, A. E. ; et al. : Experimental Investigations of Influence of Boundary 

Layers upon the Performance of MHD Generator Within a Wide Range of Hall 
Parameters. Fifth Int. Conf. on Magnetohydrodynamic Power Generation, 
vol. I, Int. Atomic Energy Agency, 1971, pp. 335-350. 

30. Nichols, Lester D. ; and Sovie, Ronald J. : Hall Current Effects in the Lewis MHD 

Generator. NASA TM X-2606, 1972. 

31. Keuhne, W. D. ; and Kolb, G. : Angewandte Magnetohydrodynamik-Heft 9. 

Berechnung der Kenngrossen von Arbeitsgasen des MHD-Verbrennungsgas- 
Generators (Applied Magnetohydrodynamics. Vol. 9 - Calculation of the Char- 
acteristic Quantity of the Working Gas of the MHD -Combustion Gas Generator). 
JUL-874-TP, Institut fur Technische Physik (Julich, West Germany), 1972. 

32. Gaponov, I. M. ; et al. : Electric Conductivity of a Plasma Composed of Natural Gas 

Combustion Products Seeded With Potash. High Temp. , vol. 10, no. 5, 1973, 
pp. 992-993. 

33. Rosa, Richard J. : Magnetohydrodynamic Energy Conversion. McGraw-Hill Book 

Co. , Inc. , 1968. 

34. Raeder, J. ; and Bu'nde, R. : Basic Equations and Fundamental Data for Combustion 

MHD Generators. Max- Planck IPP-IV-60, Institute for Plasmaphysics, Garching 
(West Germany), 1973. 


25 



TABLE I. - EXPERIMENTAL DATA FOR CALCULATION 


OF GENERATOR COEFFICIENT 


Generator type 

Refer - 

Hall 

Hall 

Gas flow 

Magnetic 

Tan a a , 

Generator 


ence 

param- 

electric 

velocity 

field 


<p 

coefficient 




eter, 

field, 

u 

strength, 






/3 

- E x 

M/sec 

B 







volt/M 


Tesla 



Combustion 

12 

1.45 

772 

1630 

1. 9 

0 

0. 172 


13 

1. 10 

443 

641 

2.5 



.246 


14 

4. 15 

3360 

2250 

3. 1 



. 116 


15 

2. 8 

(b) 

2250 

1.9 



b . 271 


16 

1.45 

1480 

1630 

1.9 


578 

.319 


17 

2.7 

1580 

1000 

2.2 

1. 0 

.39 


18 

1.3 

4000 

2500 

1.9 

1 

. 0 

.725 


19 

2.4 

1210 

1000 

3.4 

0 

. 148 


20 

1.7 

1250 

1400 

2.4 

0 

c . 38 


21 

1.9 

(d) 

720 

3. 3 

0 

d . 43 

Inert gas 

22 

2.9 

4560 

1430 

3.7 

0 

0. 297 


23 

24 

4.0 

4300 

4500 

1. 13 



.211 


25 

1. 1 

700 

1050 

. 65 



.93 




1. 15 

940 

1030 

1.03 



.77 




2.0 

1750 

1030 

1.78 



.47 




2.4 

2430 

1030 

2.72 



. 36 




2.3 

1780 

1050 

1.35 



. 54 




3. 85 

2430 

1050 

2.26 



.28 




5.0 

2540 

1050 

2.94 



. 17 


26 

5. 1 

1030 

3600 

1.4 



.04 


27 

3.9 

620 

985 

2.2 



.073 


28 

4. 5 

750 

1140 

1.4 



. 105 


29 

2.8 

— 

— 

— 



.234 




5.0 

— 

— 

— 



. 140 


30 

2.3 

45 

215 

. 19 



.472 




4.2 

88 

215 

. 34 

T 

.289 


a See fig. 2 for definition of a. 
b K H (0) = — 1 . 


K r = 


Uk h = 


- E x 

(uB + E y )’ 
(l + d 2 )P 


- 1 


I ctu 2 B 2 
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TABLE H. - MOLE FRACTIONS FOR STOICHIOMETRIC 
HYDROGEN AND OXYGEN COMBUSTION 
WITH CESIUM SEED 


[Seed: 4% CsOH/2 moles H 2 0. ] 


Species 3, 

Chamber 

Throat 

Exit 

Exit 

Exit 



o 

Pressure, N/m 



1. 034X10 6 

6. 006X10 5 

2. 811X1 0 5 

1. 39 6X1 0 4 

6. 370X10 4 



Temperature, 

K 



3320 

3178 

2992 

2832 

2663 



Mach number 



0 

1. 0000 

1. 5820 

1.9993 

2.4085 


Mole fractions 

Cs 

0. 00322 

0. 00344 

0. 00373 

0.00397 

0. 00419 

Cs + 

. 00032 

.00030 

.00027 

.00023 

. 00019 

CsO 

.00023 

.00019 

.00014 

.00011 

. 00007 

CsOH 

. 01327 

.01337 

.01350 

.01363 

.01381 

E 

. 00024 

.00023 

.00022 

.00020 

. 00018 

H 

.04222 

.03671 

. 02946 

.02330 

. 01706 

ho 2 

. 00006 

.00004 

. 00002 

. 00001 

. 00001 

H 2 

. 12461 

. 11562 

.10259 

.09017 

.07 587 

h 2 o 

. 66213 

.69146 

.73199 

.76861 

. 80849 

H 2°2 

. 00001 

.00001 

0 

0 

0 

O 

.01911 

.01617 

. 01244 

.00943 

.00655 

O' 

. 00001 

. 00001 

0 

0 

0 

OH 

.09388 

.08353 

.06961 

. 05747 

.04477 

OH" 

. 00008 

. 00006 

. 00004 

. 00003 

. 00002 

°2 

.04059 

.03885 

.03597 

.03283 

.02 880 


Additional species considered but whose mole fractions were less 
than 5X10~® Cs(S), Cs(L), CsOH(S), CsOH(L), CsOH + , Cs 2 , 
cs 2 o, Cs 2 o 2 h 2 , h + , h', h 2 o(s), h 2 o(l), o + , oh + , o 2 - , o 3 . 
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Figure 3. - Best apparent hall parameter data. 
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Figure 4. - Cross sections for momentum transfer. 




Figure 5. - Comparison of water cross sections used by various 
researchers in MHD field. 
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Oxygen to hydrogen pressure, 
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(a) Area ratio, 1.3; combustion pressure, 
1. 034 megapascals. 


2 4 6 8 

Seed fraction, mole 

(b) Area ratio, 1.3; oxygen to 
hydrogen weight ratio, 8. 


2 4 6 8 


(c) Oxygen to hydrogen weight 
ratio, 8; combustion pres- 
sure, L 034 megapascals. 


Figure 6. - Influence of percent seed on optimum values of conductivity and power densities. Magnetic field strength 
(applicable to P e plots only; see text), 10 tesla. 
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(a; Combustion pressure, 0.5171 megapascals. (b) Combustion pressure, L034 megapascals. (c) Combustion pressure, 4. 137 megapascals. 

Figure 7. - Influence of oxygen to hydrogen weight ratio on conductivity and power densities. Seed fraction, 3. 5 mole percent; magnetic field 
strength (applicable to P £ plots only), 10 teslas. 
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Figure 8. - Influence of Mach number on power densities. Seed fraction, 3. 5 mole percent; oxygen to 
hydrogen weight ratio, 7.5; magnetic field strength (P e plot only), 10 teslas. 
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Figure 9. - Combustion pressure for maximum effective power 
density versus magnetic field. Oxygen to hydrogen weight 
ratio, 7. 5 ; seed fraction, 3. 5 mole percent. 
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Figure 10. - Mach number for maximum effective power density 
versus magnetic field. Oxygen to hydrogen weight ratio, 7. 5; 
seed fraction, 3. 5 mole percent. 
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Figure 11. - Inlet static pressure for maximum effective power 
density versus magnetic field. Oxygen to hydrogen weight 
ratio, 7.5; seed fraction, 3.5 mole percent. 
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